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HOMOGENEOUS TURBULENCE IN A COMPRESSIBLE FLOW 
by Robert G. Deissler 
Lewis Research Center 


SUMMARY 

Two-point correlation equations for turbulence with a uniform strain rate 
in one direction are constructed by starting with the compressible Navier- 
Stokes and continuity equations. Compressibility is retained in the mean flow, 
but the problem is simplified by assuming that the turbulent eddies behave in- 
compressibly . In order to close the system of equations, the weak turbulence 
approximation (triple correlations neglected) is used. The equations are con- 
verted to spectral form by taking their Fourier transforms. Solutions are ob- 
tained for the turbulent velocity variances in the longitudinal and transverse 
directions and for the spectra of those quantities . The transverse turbulent 
intensity (square root of variance over local mean velocity) varies in a com- 
plex manner that seems to approximate experimental heat transfer results in 
highly heated tubes. 

A transfer term associated with the mean strain rate is found to transfer 
energy to either higher or lower wave numbers, depending on whether the strain 
rate is negative or positive. Because of this energy transfer the peaks of the 
energy and dissipation spectra are widely separated for large negative rates of 
strain, even though triple correlations have been neglected; however, there is 
still considerable overlap of the two spectra. 


INTRODUCTION 

One problem of considerable interest in the theory of turbulent flow is 
the effect of a mean strain on the turbulence. Since the strain is accompanied 
by a velocity change in the flow direction, this effect is sometimes called an 
effect of acceleration. The case of turbulence with a suddenly applied mean 
strain, in which viscous and inertia effects produced by the turbulence are 
negligible, was considered, for instance, by Prandtl (ref. l), Taylor (ref. 2), 
Ribner and Tucker (ref. 3), and Batchelor and Proudman (ref. 4). A different 
approach wherein viscous effects are retained and inertia effects are neglected 
only insofar as they affect triple correlations was used by Pearson (ref. 5). 
Two-point correlation equations in which uniform mean strain rates are in- 
cluded were obtained from the incompressible Navier-Stokes and continuity equa- 
tions . 
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In the preceding analyses, a strain in the flow direction was accompanied 
by transverse strain, as in the distortion produced by a sudden stream con- 
traction. The problem analyzed herein, on the other hand, considers the tur- 
bulence in a stream of uniform cross section, in which a uniform rate of strain 
is produced in the flow direction by density changes. This type of strain 
might be of importance in connection with heat transfer to or from a gas stream 
at high fluxes or in a flow with heat liberating or absorbing chemical reac- 
tions. The model might be closely approximated in an experiment in which a 
number of very small heated or cooled tubes run longitudinally in the stream. 
The turbulence would be initially produced by flow of the stream through a 
grid, and the effect of the strain rate on the decay of the turbulence would 
be studied. 

The analysis starts from the compressible Navier-Stokes and continuity 
equations. Erom those equations two-point correlation equations for turbulence 
with a uniform mean rate of strain in one direction are obtained. Compressi- 
bility must, of course, be retained in the mean flow, but density fluctuations 
are neglected and the turbulent eddies are assumed to obey the incompressible 
continuity equation. The required correlation equations will be obtained in 
the following section. 


BASIC EQUATIONS 

The compressible Navier-Stokes and continuity equations for a constant- 
viscosity gas are usually written as follows: 



where the subscripts can take the values 1, 2, or 3. A repeated subscript in 
a term indicates a summation of terms, with the subscripts successively taking 
on the values 1, 2, and 3. The quantities up and u^ are instantaneous 
velocity components, xp is a space^coordinate, t is the time, p is the 
density, |i is the viscosity, and p is the instantaneous pressure. (All 
symbols are defined in the appendix. ) 

Equations (l) and ( 2 ) give 

It (p ^ } + (pSk5i) 

The instantaneous velocities and pressure in equations (3) and ( 2 ) can be 
broken into mean and fluctuating components. Set up = Up + up, and 
p = P + p, where the capitalized symbols refer to mean quantities, and the 
lower case symbols refer to fluctuations. Density fluctuations are neglected. 


3 oxi 


dXi 
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If averages are taken and the averaged equations are subtracted from the un- 
averaged ones. 


sl + (P u i u k) + ^ P u k u i ) + ^ (pu ± u k ) - ^ (pu iUk ) 



Since continuity is satisfied for the mean as well as for the unaveraged flow, 
equation (2) becomes 


St + W - 0 (5) 

Then, from equations (2) and (5), 

S^ < pu k) a 0 (6) 

Also, since fluctuations are assumed to behave incompressibly, the incompressi- 
ble continuity equation is used: 


Su k 


0 


( 7 ) 


Comparison of this equation with equation (6) indicates that its use implies 
that changes in mean density are small over the distance for which a velocity 
fluctuation occurs. This distance can also be considered as the mixing length 
or scale of the turbulence or the distance over which velocities are appre- 
ciably correlated. Equation (7) allows dropping of the last term in equa- 
tion (l). Use of equations (5) and (6) then gives 


Su ± Sup SCpu^) S(pupu k ) S 2 up 

P Tt' + pu k S^ + pU k S^ + = "^p + P 


( 8 ) 


Equation (8) applies at a point P in the fluid. At another point P' sepa- 
rated from P by the vector r (see fig. l). 


P' 


Sul SU Sul S( p ' uA u t ) 

■St + p,u k 5^ + P' u k sf + 


a(p ' u ,H ) v . a2u ,i 

Sx' Sx' p Sx'Sx' 

ox k ox j ox k ox k 


(9) 


If equation (8) is multiplied by p 'uj and equation (9) by pup, the 
two equations are added, and average values are taken. 
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du.u! 3U. dU. du.ul 

PP' -gp- + PP’ ^ + pp' u ± u£ ^ + pp'U k 


, ^ u i u i ^(pp'UiU^Ulj) ^(pp'UiU^) 

+ «>’ u k -sf + — * — LJ - 


5x k 


"5x! 


, dpui duj_p' , , d^u-jui L dniju* 

= -P' - P - + P V s \ J + PH s..Vs ry 




dx k dx k dx k dx k 
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( 10 ) 


u' 



Figure 1. - Vector configuration for two-point 
correlation equations. 


where the fact that quantities at one point are 
independent of the position of the other point has 
been used. As mentioned previously, use of the 
incompressible continuity equation for the fluctu- 
ations (eq. (7)) implies that changes in mean den- 
sity are small over the distance for which ve- 
locities are correlated. Therefore, p* is re- 
placed by p, and in that way a considerable 
simplification in equation (10) is obtained. Also, 
the turbulence is assumed to be weak enough for 
triple correlations to be neglected. Equation 
(10) then becomes, on introducing the new vari- 
ables r k = x£ - x k and (x k ) m = (l/2)(x k + x£) 
(see fig. l). 


Su n -u!s dU,- BU- Ut, + IXu chi-iu- . du-tu- 

u k u. UiUi4 + ^^ + - U k ) -jiJ. 


m 


_1 

2p 


dpu 


m 


du.p ! 

tkr 

m 


x /dpui 

+ p 1 st - "5?7 



, H , 2p ^ u i u ,i 

2 P ^(x k ) m S(x k ) m P ^r k dr k 


( 11 ) 


where the following transformations were used: 

S 1 d S 

5^-2 snjr- 5^ 

m 

d _ 1 5 a 

^ = 2 bl^) ^ 
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For the problem considered herein, the velocity gradient (strain rate) in the 
xq-direction dUq/dxq is uniform, and the other velocity gradients are zero. 
Also, it is shown in reference 6 that the effects of inhomogeneities on the 
decay of the turbulence downstream of a grid are, for cases of practical in- 
terest, usually negligible. Equation (ll) then becomes, for steady state at a 
given distance downstream from the turbulence generator, 



where (Up + Up)/2 was replaced by (Up) , the value of Up at (xq) , and 6p< 
is the Kronecker delta . m m 


The pressure-velocity correlations are next considered. Taking the 
divergence of equation (4) and using equations (6) and (7) give 



Since the mean flow is in the xq-direction and changes only in the xq- 
direction, the first term on the right side of equation (13) becomes, with the 
use of equations (5) and (7) and the assumption of steady state. 





= 0 


Multiplying equation (13) by uj, averaging, and introducing the variables 
(xq) and rq then give, for weak turbulence. 



Equation (14), h owev er , is the same as equation (12) in reference 6. As in 

that reference pul is assumed to be zero at an initial (x-, ) . This assump- 

J - 1 - m 

tion is consistent with an assumption to be made l ater that the turbulence is 

initially isotropic . In the preceding reference puj is shown to be zero 

throughout the field if it is zero at an initial (xq) : that is, the dynamics 
of the turbulence will not cause nonzero values of puj to arise. Similarly, 
there will be no inconsistencies in also taking uqp' to be zero. Equation 
(12) then becomes 
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In order to simplify the notation, let 


(Up) = U and (xp) = x 
•*- m x m 


By continuity, 


PU = Pg U g 


or 


Also, 


P U 

n - lS-£ 

p - U 


U = u g + 



(16) 


(17) 


where the subscript g designates the plane where x = 0 which, as usual, is 
taken as the plane of the grid. If 


and 



dU 


u e 

X = — + 
s 


x 


v 


= d 


equation (15) becomes 


5 ii u i u j 


5 . , u . u' + X 
jl i 1 


^ u i u j 

dX 


du i u j 2v w 

ri ^T = \ x ^k 


u i u i 


(18) 

(19) 


( 20 ) 


Equation (20) is the correlation equation for u-u-. It can be converted to 

spectral form by introducing the three-dimensional Fourier transform cpp. de- 
fined by 


( 21 ) 


# 


(Pile** 1 dx 


^-f 

•'-03 

where dx = dx-^ dx 2 dXj. Then, 

Taking the Fourier transform of equation (20) results in 
^ii U g U r U s / ^ii ' 

U g “3x = " 8 il f "Pij - 6 jl # ^il + # K 1 + ^ij 


2vK 2 cp i j 


( 22 ) 


(23) 


If X/U g is interpreted as a time, the terms on the right side of equa- 
tion (23) give the contributions of various processes to the time rate of 
change of j • The last term is, of course, the usual dissipation term. The 
next to the last term can be interpreted as a transfer term, for if r = 0 in 
equation (22), the quantity (k^ dcp-y/ck-j_) + Tij* when integrated over wave 
number space, gives zero contribution to the rate of change of UqUj . It can, 
however, transfer energy between wave numbers. This process is similar to that 
previously discussed in connection with shear-flow turbulence (ref. 6). The 
remainder of the terms on the right side of equation (23) can evidently be 
interpreted as production terms. Both the production and transfer terms drop 
out of the equation when there is no strain, since X goes to infinity when 
s goes to zero (eq. (19)). 


SOLUTION OF SPECTRAL EQUATION 

Equation (23) is a first-order equation and can be solved by available 
methods (see ref. 7). For cpqq, the component of cpqj in the direction of 
mean strain. 



(24) 


where f(fCpX) is a function of integration that depends on initial conditions. 
In order to evaluate f it is assumed that the turbulence is isotropic at 
X = Xq (but not at other values of X). This assumption implies that, for 
weak turbulence. 


(<P,J = 


J 0 


( K% ij - Vj) 


(25) 


^'O 12k 2 

where Jq is a constant (see ref. 6, eq. (43)). This equation is in agreement 
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with the work of Batchelor and Proudman (ref. 8). Evaluation of f in equa- 
tion (24) by substituting equation (25) at X = X 0 gives 

r _ 

. .... \2 

„2 , , 2 . ' 


f(KiX 0 ) = x 0 -2- (k| + K§)exp«{ 

12it^ 1 


2vX r 


U 


g 


(*1 X q ) ; 

•v-2 

L X 0 


- K? 



or 


f(k 1 X) = — x 0 (k| + K 3 )exp 


12it 


2vX Q /k|x 2 

TAT 


K 2 - K 2 
k 2 k 5, 


Then 
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- - >8>»4c ■ ■.>[* * •©’ ijr^ 


^ 12it 2 ^g + sx 


"g 


(26) 


where the definition of x (eq. (19)) was used. Similarly, cp , a component 
normal to the direction of mean strain, is 2 


0 

^ 22 “ ~~2 
12jc^ 


u g + sx h + 2/ u g + sx \ 2 ] 

Ug + -0 L 3 % + 3 V J 


X exp^- — K 2 (x - X 0 ) 
g 


i + s (hf x - x ° 

\K) Ug + SX C 


(27) 


Because of axial symmetry, cp 33 is of the same form as cp 22 and need not be 
considered separately. 

In order to integrate over wave number space, it is convenient to intro- 
duce spherical coordinates as follows : 


K-^ = k cos 6 


k 2 = K cos cp sin 6 


= K sin cp sin 6 


Equation (2l) then becomes, for r = 0, 

“T 1 ! - 1 *u 


(28) 


(29) 
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where 


t 


,-// 


2rt 

, . cp. ,K 2 sin Q dcp d0 

ij J J ij Y 

0 0 


(30) 


The quantity ^ is the spectrum tensor as defined by Batchelor (ref. 9). 

In calculating u 2 and u 2 it is convenient to carry out the integration 

with respect to cp, K, and 6 in that order. Then, from equations (26) 
to (30), 


2 J 0 


Ul 48-v/2it" \x - x 0 J Vu + sx 


Ug/v 


V 5 /^ P + sxqV V 2 


g 


(31) 


and 


2 2 J 0 


u 


u 2 = u 3 = 


,/v \ 5 / 2 [s ( x + x 0 ) + 2U g ](Ug + sx) 1 / 2 


48V2rt 


x - Xx 


& r _ 

2(U g + sx 0 ) 3 / 2 


(32) 


For the spectra of the velocity variances there is, in dimensionless form. 


.* b K*e -2K: 

tn = 


*2 


11 12ita 


[W!LlJ A erf (K* -\fzai) + K*e- 2aK * 2 ] 


(33) 


and 


■*2 


* 4 * _ kV 2ic 

22 ' *53 2toab 3 


jy2na erf(K* -J2a) + 2k* 2 J 


+ (b 2 - 2)k%“ 2aK 


*2*^ 


for a > 0 and 




* _ b K 

11 “ 


+ F( y^a k*) ♦ k*] 


(34) 


(35) 


and 


»* = ** - «*exp[- 2K *2(l + a)] pttex* 2 - 1) + 2 + (b 2 . 2)k *] 

53 24:tab 5 L V"2a J 


( 36 ) 


for a < 0, where 
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and 



F(co) 



(42) 


Values of F(<o) are tabulated, for instance, by Miller and Gordon (ref. 10). 


The transfer term in the spectral equation (the next to the last term in 
eq. (23)) can be integrated over all directions in wave number space by re- 
placing cp^j in equation (30) by that term. The longitudinal component of the 

integrated transfer term is, then, in dimensionless form, 


rn* _ 

i ll" 


bK*e~ 2K 


*2 


12na 
for a > 0 and 
,* 


4aK* 2 - 5 
4a 


^2ita erf(K*-^/2a) + 3 K*e“ 2aX 


*2 


m* _ _ bK exp[-2K 2 (l + a)] 
12na 


4aK* 2 - 5 
a/ - 2a 


F( -y^2a K*) + 3 K 


+ abi|r*^ 


+ abi^h 


(43) 


(44) 


for a < 0, where 
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(45) 


* 


, 2 (x 


* ■ *o )3s 

T H ~ VrrZfrr" . ~ T 11 


J 0 U|(U g + sx) 


A discussion and plotting of the results for the various spectra and velocity- 
variances will be given in the following section. 


DISCUSSION 
Velocity Variances 

— ry — ry — rr 

The velocity variances Ug, and are given by equations (3l) and 

(32). For zero strain (s = 0) these expressions reduce to the usual expression 
for isotropic turbulence in the final period, if (x - Xg)/Ug = t - tQ. In 
order to plot equations (3l) and (32) for s / 0, the sign of xq must be 
known. It appears that xo must always be negative, since equations (3l) and 
(32) indicate that the point x = xq represents a virtual origin for the tur- 
bulence, where the energy is infinite. Since the grid at x = 0 can produce 
only a finite turbulent energy, the virtual origin must lie upstream of the 
grid, or x Q must be negative. 

Dimensionless plots of equations (3l) and (32) are presented in figure 2 
for two values of x*. The velocity variances at a given x with strain are 

divided by their values at the same position without strain, with the mean 

velocity set at Ug. For positive strain rates the longitudinal component u| 
(in the direction of strain) is reduced below the value it would have for no 
strain, whereas the transverse components are increased. These results are 
qualitatively similar to those obtained for a sudden contraction in an incom- 
pressible flow (see ref. 9). For negative strain rates the opposite trends 
occur; that is, the longitudinal component is increased by the strain, whereas 

the transverse components are decreased. The total energy UjU^/ 2 = (u^ + 2u|)/2 

is increased by both positive and negative strain rates; that is, positive and 
negative strain rates both feed net energy into the turbulent field but through 
different components . 

As s x approaches -1 (or sx -* -Ug), and thus Uj_Uj_ grow in- 

definitely large (see eq. (3l)). The symbols s* and x* are defined in 
figure 2 . The point sx = -Ug represents the condition where the fluid is 
compressed into zero volume and is, of course, unattainable. 

For values of sx somewhat larger than -Ug, however, the magnitudes of 

u 2 and u i u i still increase rather than decay as x increases, or the energy 

fed into the turbulent field by the negative straining action is greater than 
that dissipated. It can be shown from equation (3l) that the region in which 

u 2 increases with x is given by -l/x* < s* < -(5/8)/[x* + (3/8)]. 

~ 2 " / 2 

Calculating the square of the transverse turbulent intensity u 2 /U , 
where U is the local mean velocity, is of some interest. In a sudden con- 
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Shortening 


Elongating 


2 

2 

u 2 

U 3 

2) 

(72) 

2L 

l U 3/n 


CM 

=J 

(a) x* ■ x/(-xq) = 1. 

"n, 10 1— 

— •? r 

CM CM I— 

[ | (Shortening I Elongating 

1 1 1 


S* = S(-X Q )/U g 
(b) x* = x/(-Xq) = 3. 

Figure 2 . - Ratio of velocity variances with strain to those without strain as 
function of dimensionless mean strain rate. 


traction, results for an incompressible flow indicate that, although Ug in 

creases, u 2 /U 2 decreases through the contraction because of the increase in 
U (ref. 9). In the case presented herein it is desired to calculate 

^u| at x with strain^(U^ at x with strain) 

Ul at x without strain and with U = 
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where equation ( 32 ) and the relations 
U = U 0 + s(x - xq) and U g = U 0 - sxq 
were used. 


A plot of equation ( 46 ) is given 
in figure 3. For small values of x* 

(ng. 3(a», 

tends to be greater than 1 for positive 
strain; that is, the positive strain 
pate tends to increase the local trans- 
verse intensity of the turbulence. For 
larger values of x*, the ratio is pre- 
dominantly less than 1 except for large 
values of s*. As x* approaches in- 
finity, equation (46) shows that the 
ratio of local transverse intensity 
with strain to that without strain goes 
to zero. 

Comparing these results for posi- 
tive strain rates with those for expand- 
ing flow in a highly heated tube may be 
worthwhile, although the heated- tube 
case is more complicated because of 
radial variations of velocities and 
properties . Experiments by Weiland 
(ref. 11; see fig. 7) and Taylor 
(ref. 12) indicate that the heat trans- 
fer in an expanding flow can be con- 
siderably different than might be ex- 
pected for a nonexpanding flow. In the upstream portion of the tubes the heat- 
transfer coefficient tended to be higher than calculated for a normal entrance 
region. In the middle portion the heat-transfer coefficient dipped below the 
normal value, and near the exit of the longer tube it again rose. The experi- 
mental curves are qualitatively similar to those in figure 3(a), except that 
the latter do not rise near the exit. This difference can probably be attri- 
buted to the fact that the strain rate was not uniform along the length of the 
tube in the experiment. The dimensionless heat-transfer coefficient or Stanton 
number is closely related to the transverse turbulent intensity. In fact, some 
order of magnitude arguments based on the momentum - heat -transfer analogy can 
be used to show that the ordinate in figure 3 should be approximately propor- 
tional to the square of the ratio of local Stanton number with strain to that 
without strain. The results presented herein thus seem to be in qualitative 
agreement with the experimental heat-transfer results. Although, as mentioned 
previously, the heated-tube case is more complicated than the case considered 
herein, the results appear to show, at least, that appreciable effects of axial 
position on heat transfer at distances far downstream from the entrance might 
occur. 

For negative strain rates (fig. 3(b)), the trends are, in general, oppo- 
site to those for positive strain rates; for small values of x*. 



Figure 3. - Square of ratio of local transverse inten- 
sity of turbulence with strain to that without strain 
as function of dimensionless axial position and di- 
mensionless strain rate. 
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Figure 4. - Dimensionless longitudinal transfer spectra associated with 
mean strain at dimensionless axial position x" of 1. (Arrows point 
to ordinate scales to be used for various curves. ) 




is less than 1, and for larger values of x*, it becomes 


greater than 1. If the momentum - heat-transfer analogy is applied as it was 
for positive strain rates, it is seen that negative strain rates could either 
decrease or increase the heat-transfer coefficient. Negative strain rates 
occur, for instance, in the gas flowing through a cooled tube. 


Turbulence Spectra 

The dimensionless transfer spectra are next considered. The transfer 
term in the spectral equation for the longitudinal component of energy is 
plotted in figure 4 from equations (43) and (44). These curves and the suc- 
ceeding ones are for a value of x* of 1. Curves for other values of x* are 
qualitatively similar . For negative strain rates the transfer term is negative 
for low values of K* and positive for large values, or the energy transfer is 
from small to large wave numbers. The total area under each curve as K* goes 
from zero to infinity is, of course, zero. The curves are similar to those 
obtained for the energy transfer associated with triple correlations at higher 
Reynolds numbers (ref. 13). In the case considered herein, where the triple 
correlations are neglected, the energy transfer is caused by the mean strain. 
For positive strain rates the energy transfer is in the opposite direction, 
that is, from large to small wave numbers. The energy transfer can thus be in 
either direction, depending on whether the strain rate is positive or negative. 
The transfer spectrum associated with the transverse energy, that is, the plot 
of Tg 2 , although not shown, is similar to that of Tqq, and the direction of 
the energy transfer corresponds to that for T^. 

-¥r 

Dimensionless spectra of the velocity variances, or values of tn and 
1 ^ 22 , are plotted in figures 5 and 6. For negative strain rates (figs. 5(a) 
and 6(a)), the shapes of the spectra are changed considerably because of the 
transfer of energy into the high wave number region as discussed in the pre- 
ceding paragraph. The excitation of the high wave number portions of the 
spectra causes them to become asymmetric, as compared with the almost symmetric 
spectrum for no strain (s* = 0). In the case of the curves for negative 
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Figure 5. - Dimensionless longitudinal energy spectra at dimensionless 
axial position x* of 1. 


. o 

strain rates overlap because, although contributions to u^ with strain are 

generally lower than those without strain, the opposite trend occurs at high 
wave numbers because of the energy transfer into that region. For positive 
strain rates (figs. 5(b) and 6(b)), where the energy transfer is from large to 
small wave numbers, the shapes of the spectra are affected to a lesser extent by 
the energy transfer. Apparently, the effect of the energy transfer for positive 
strain rates is to move the peaks of the curves to the left, particularly those 
for f 22 . 

Figure 7 is a summary plot for energy, dissipation, and transfer spectra 
for a high negative strain rate (s* = -0.99). The peaks of the energy and 
dissipation spectra are quite widely separated and occur at values of K* of 
1.3 and 10, respectively. This separation is similar to the separation that 
occurs at high Reynolds numbers without mean strain (ref. 13) and is a con- 
sequence of the asymmetrical shape of the energy spectrum. The entire energy 
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and dissipation regions, however, cannot be said to be widely separated, as is 
required for Kolmogorof f ' s hypothesis to be valid (see ref. 9) because there is 
considerable overlap of the spectra. In an attempt to obtain less overlap, a 
still higher value of negative strain rate (s* = -0.999) was tried. Although 
the peak of the dissipation spectrum moved to about 32, while that of the 
energy spectrum was unchanged, the overlap of the regions was not decreased 
because ordinate values for the energy spectrum in the high wave number region 
were increased. It seems that a nearly complete separation of the energy and 
dissipation regions cannot be obtained with the present model. 


Lewis Research Center, 

National Aeronautics and Space Administration, 
Cleveland, Ohio, February 17, 1965. 
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APPENDIX - SYMBOLS 


a 

b 

F 

J 0 

P 

P 

r 

s 


defined by eq. (37) 
defined by eq. (38) 
defined by eq. (42) 

constant that depends on initial conditions 
mean pressure; also, a point in the fluid 
pressure fluctuation 
defined in fig. 1 

^ = au 


longitudinal transfer term 
t time 

U,(Ui) m mean longitudinal velocity component 

u velocity 


X 


X;(xi) m 

8 ij 

K 

M- 



X 


s 

longitudinal space coordinate 

Krone cker delta 

angular coordinates (see eq. 

■wave number 

viscosity 


(28)) 


v p/pg 

P density 

three-dimensional Fourier transform of u-^uj defined by eq. (2l) 
\|/ - i spectrum tensor defined by eq. (30) 

Subscripts : 


g at plane where x = 0 (plane of grid) 

m at point P m (see fig. l) 
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0 


at virtual origin of turbulence where turbulent energy would be in- 
finite (It is assumed that turbulence is isotropic at xq and that 
strain begins to act there . ) 

1 in longitudinal direction 

2,3 in transverse direction 

Superscripts : 

~ instantaneous 

' at point p' (see fig. l) 

* on dimensionless quantities defined by eqs . (39), (40), and (45) 
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